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Abstract 

We consider "Orientifold field theories", namely SU(N) gauge the- 
ories with Dirac fermions in the two-index representation at high tem- 
perature. When iV is even these theories exhibit a spontaneously bro- 
ken Z2 centre symmetry. We study aspects of the domain wall that 
interpolates between the two vacua of the theory. In particular we 
calculate its tension to two-loop order. We compare its tension to 
the corresponding domain wall in a SU(N) gauge theory with adjoint 
fermions and find an agreement at large- N, as expected from planar 
equivalence between the two theories. Moreover, we provide a non- 
perturbative proof for the coincidence of the tensions at large- N. We 
also discuss the vacuum structure of the theory when the fermion is 
given a large mass and argue that there exist N — 2 metastable vacua. 
We calculate the lifetime of those vacua in the thin wall approximation. 



1 Introduction 



When a global discrete symmetry is spontaneously broken there exist 
domains walls that interpolate between the various degenerate vacua. An 
important example is the case of the deconfining phase of gauge theories 
that admit a non-trivial centre symmetry. Due to asymptotic freedom, the 
high temperature phase is controlled by a weak gauge coupling, hence it is 
possible to calculate the domain walls tension order by order in perturbation 
theory. The tension of the pure SU (N) Yang-Mills theory was calculated at 
the one- loop order in [1, 2]. The calculation was extended to two-loop and 
three-loop order in [3] and [4]. Up to two loop order the tension exhibits 
a Casimir scaling. A deviation from Casimir scaling was found at three- 
loop order [4]. Following [5] the tension of the deconfining domain walls in 
Af = 4 SYM was calculated at both one-loop and strong coupling in [6] and 
extended to two loop order in [7]. 

In this paper we focus on "orientifold field theories" (or OrientiQCD), 
namely on SU(N) gauge theories with Dirac fermions in either the two- 
index symmetric or the two-index antisymmetric representations. When N 
is even the centre symmetry is Z 2 [8]. Therefore there exists a domain wall 
that interpolate between the two vacua of the theory. The wall tension is 
0{N 2 ). 

The main motivation for our investigation is the large- N equivalence be- 
tween SU(N) theory with Nf Dirac fermions in the two-index (either sym- 
metric or antisymmetric) representation and SU(N) theory with Nf Majo- 
rana fermions in the adjoint representation [9, 10]. Planar equivalence states 
that OrientiQCD and AdjQCD (adjoint QCD) become equivalent in a com- 
mon sector of bosonic charge conjugation (C-parity) invariant states if and 
only if C-parity is unbroken in OrientiQCD [11]. To present, both analytic 
and numerical works focused on light states, namely on glueballs and mesons 
as well as the quark condensate. Heavy objects, such as domain walls, were 
not yet investigated. 

In the deconfining phase the equivalence between the two theories holds 
in the common vacua of the the theories, namely when the Polyakov loop 
expectation value is either 1 or —1 [11]. Apriori it is not clear whether the 
equivalence holds for the domain wall, since it interpolates between different 
vacua. In addition the domain wall tension is 0(N 2 ), namely it is a "heavy" 
object, and it is not obvious that large- N equivalence should hold for such 
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heavy objects. 1 

We carry out a two-loop calculation of the domain wall tension for an "ori- 
entifold theory" with Nf massive fermions. The result is that the tension of 
the wall matches the tension of the corresponding wall of the adjoint fermions 
theory in the large- N limit, namely that the leading 0(N 2 ) contribution is 
the same in both theories. 

In addition to the explicit two-loop calculation we provide a non-perturbative 
proof for the large-iV coincidence of the tensions. The proof is based on 
the orientifold planar equivalence as formulated in terms of the coherent 
states [12, 11]. The main ingredient will be the reformulation of the domain 
wall tension at finite temperature as the tension of a Dirac magnetic string 
at zero temperature but on a compactified space. 

Another part of our paper concerns metastable vacua in "orientifold" the- 
ories. When the two-index fermions are given infinite mass they decouple 
and we end-up with the Yang-Mills deconfining vacua, namely N degenerate 
vacua which correspond to the spontaneously broken Z N centre symmetry of 
the pure Yang-Mills theory. Therefore at large (but finite) value of fermion 
mass there should be two real vacua and N — 2 metastable vacua. We con- 
firm this scenario by an explicit one-loop calculation of the potential for the 
Polyakov loop. We then discuss the lifetime of the metastable vacua in the 
thin wall approximation [13, 14]. 

The paper is organized as follows: in section 2 we describe the general 
framework of the paper. In section 3 we present the results of the two-loop 
wall tension calculation (leaving the technical details to the appendix). In 
section 4 we provide the non-perturbative proof for the large- N coincidence 
of the wall tensions. Section 7 is devoted to the metastable vacua and their 
lifetime. In section 6 we discuss our results. 

2 Two guises for the 't Hooft loop 

This section provides a platform for the perturbative approach in section 

3 and the non-perturbative proof in section 4. We will review the original 
formulation by 't Hooft [15] for a spatial loop in a theory with non-zero tem- 
perature (so fermion fields will be antiperiodic in the Euclidean time direc- 
tion), in terms of a discontinuous gauge transformation. This one is strictly 

1 0{N 2 ) domain walls are described by wrapped NS5 branes in the dual string theory 

[6]. 
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related to the usual formulation as an electric flux loop operator, which will 
be useful for the perturbative proof (to low orders) of the equivalence. 

Rotating one of the loop sides in the temporal direction is the starting 
point for a Hamiltonian formulation, where one studies the propagation of 
a magnetic flux tube at zero temperature in a 3D space with one compact 
periodic dimension. The fermion fields are still antiperiodic in that direction. 
This Hamiltonian formulation will be useful for the non-perturbative proof 
of the equivalence in section 4. 

2.1 Continuum formulation 

We start with the Hamiltonian formulation on a 3D space (infinite in every 
direction) 2 . The Hamiltonian H is built from the electric field strengths E 
and the magnetic field strengths B. 

Apart from the SU(N) gauge field there are fermionic fields ip, repre- 
sented by Hermitean N x N matrices. We are interested only in adjoint and 
two-index fermions. The theory with adjoint fermions is invariant under the 
full center Zjy, while the theory with symmetric or antisymmetric two-index 
fermions is invariant under the Z2 subgroup if N is even. The Hamiltonian 
consists of a gauge part and a fermionic part, minimally coupling fermions to 
the gauge field through the covariant derivative D(A)ip which is dip+ig[A, ip] 
or dip + ig(Aip + ipA T ) (for adjoint and two-index fermions respectively). 

Choosing an oriented surface S subtended by the oriented loop L, the 't 
Hooft loop was originally [15] defined as the unitary operator 

V k (L) = exp jl J [E b .D(A) bc + gTr^T c i>]uj c k (x)d 3 x} , (1) 

representing a time-independent local gauge transformation exp(iT a ou k (x)) , 
that has a discontinuity of z k = exp (ik^r) through the surface S. The dot 
stands for summing over spatial degrees of freedom. This operator depends 
in general on the particular choice of the surface S. However, if z k belongs to 
the symmetry subgroup of the theory, the surface can be arbitrarily deformed 
(provided that its contour is unchanged) via a regular gauge transformation. 
Therefore in this particular case, the operator V k (L) acting on physical (i.e. 
gauge invariant) states depends only on the loop L and not on the particular 
choice of the surface S. 

2 This subsection relies on the work in ref. [16]. 
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The immediate effect of V k (L) is that any Wilson loop W(C) in the 
fundamental representation with C winding once around L in a clockwise 
direction will pick up a phase z k ('t Hooft algebra): 

V k (L)W(C)V k \L) = z k W(C) . (2) 

From standard considerations [16, 7], up to a regular gauge transforma- 
tion and provided that z k belongs to the symmetry subgroup of the theory, 
the 't Hooft loop can be rewritten as an electric flux operator: 

Vfe(L) = exp dS.TrEY^J . (3) 

The numerical N x N matrix Y k equals 




N—k times k times 



generating the center group element z k l = exp(2niY k ) . 

Deep in the deconfined phase the electric flux is due to almost free sta- 
tistically independent screened gluons, and elementary arguments show that 
its thermal average falls off exponentially with the area of the loop [3, 16]: 

(Vk(L)) = ^Tr[\4(L)e-^] * Ae'^^ , (5) 

where Z(T) = Tr e~ u l T is the partition function, and the trace is taken only 
on physical (i.e. gauge invariant) states. 

2.2 't Hooft loop as magnetic correlator 

If we write the thermal expectation value (5) in the path integral formu- 
lation, the insertion of the operator V k {L) is implemented as some twisted 
boundary conditions for the gauge field in the thermal direction: 

A fl (x,y,z,r + l/T) = A^x, y, z, r) + 2irY k 6 litX X[o,L x ] (x)x[o,L„](y)S(z) . (6) 

Since we are dealing with the Euclidean spacetime, we can now interpret y 
as the thermal axis. The 3D space is identified by the coorinates (x, z, r) 
and the direction r is a spatial compact direction. The Hamiltonian in the 
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3D space is written as TL. The boundary conditions above introduce a twist 
(localized on a string) in the compact spatial direction, only in the time lapse 

From a Hamiltonian point of view, if Sk{L x ) is the unitary operator which 
creates the twist 3 , then the thermal expectation value of the 't Hooft loop 
can be written as the zero-temperature Euclidean correlator of the twist 
operator 4 : 

(V k (L)) T = (0\S k (L x y e - nL yS k (L x )\0) . (7) 

When S k (L x ) acts on a small closed loop 7 in the z — r plane that encircles 
the string identified by r = z = and < x < L x at a fixed time, a z k factor 
is produced: 

S k (L x yW( 1 )S k (L x )=z k W( 1 ) . (8) 

Notice that this is different from the 't Hooft algebra (2), since the operator 
S k (L x ) describes a open string-like singularity. Therefore S k (L x ) can be 
interpreted as the creation operator for a Dirac string carrying a magnetic 
flux of strength z k , ending in a Dirac monopole-antimonopole pair. This 
operator causes the Polyakov loop operator 

,1/T 

Vexp(i / drA T (x,z,r)) (9) 
Jo 

as a function of z to jump with the z k factor when it crosses the location 
of the Dirac string at z — 0. This is based on the operator identity for the 
(S-transformed Polyakov loops: 

rl/T 

lim V exp(i / drAf (x, z = — 5z, r)) x 

f° 

xVexp(-i dTAf(x,z = 5z,r)) = z k , (10) 

J 1/T 

the left hand side being a closed loop around the Dirac string at z = t = 0. 



2.3 Lattice discretized Hamiltonian 

The lattice regularization is particularly useful because it smooths out 
the singularity in the the boundary conditions (6). Moreover it provides a 

This operator produces field configurations which are multivalued on the compact 
direction. 

4 For simplicity we choose the vacuum to have zero energy. 
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manifestly gauge covariant formulation. We remind that in the Hamiltonian 
formalism only the 3D space is discretized, while the time (y in the notation 
of the previous subsection) stays a continuous parameter. 

Formally we just need to replace the Hamiltonian in eq. (7) with the 
lattice version (x = (x, z, r) are the integer valued coordinates of three space): 

x y u 

where the operators Un are the Wilson loops around the smallest square 
paths on the lattice (plaquettes), and Hp is the fermionic contribution to 
the lattice Hamiltonian. We do not need to specify the form of Hp. We will 
just assume that it inherits the center symmetry subgroup from its continuum 
counterpart. 

The Hamiltonian % L is transformed by S k (L x ) into the twisted Hamilto- 
nian 

U\ = S k (L x y H L S k (L x ) (12) 

using eq. (8). The operator Sk(L x ) twists with a factor z k all the plaquettes 
that are pierced by the Dirac string: 

Hi = ^£ Tr ^ x ) 2 + i E ^Tr(l-z k U D ) 

x y twisted □ 
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+ J2 MeTr(l-Uu)+H F . (13) 



nontwisted □ 

With this definition the expectation value of the 't Hooft loop becomes: 

(V fc (L)> T = (0|e-"^|0> . (14) 

This twisted Hamiltonian will be used in section 4. Note that it is Her- 
mitean for all twists, but charge conjugation invariant only if z k — ±1. 



3 Domain wall equivalence in perturbation the- 
ory 

In this section we discuss the perturbative equivalence of the two-index 
SU (N) gauge theory (often called "orientifold field theory" or OrientiQCD) 
with M = 1 SU(N) theory (also called AdjQCD) up to terms of 0(1/N). 
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Figure 1: Monopole anti-monopole pair induced by twisting the z-r Wilson 
loops (plaquettes in lattice formulation) pierced by the Dirac string in the x 
direction. This string is propagated in the y direction. 

The two-index spinor ip^ carries two SU(N) indices % and j, running from 
1 to N. They refer to the transformation property of ip l i under an SU(N) 
element U in the fundamental representation: 

^ -> (U^U T ) lj (15) 

which leads for even iV to invariance under the subgroup Z 2 of the full centre 
group Z N . 

We start with the action of the two-index theory. Apart from the gauge 
field action it contains the minimally coupled fermion action: 

^MA)^. (16) 

Because of eq.(15) the covariant derivative acts as: 

D^A)^ = + %A^ + (17) 

instead of a commutator as behooves an adjoint fermion. The 7 matrices are 
Hermitean. 

3.1 One loop approximation to the tension 

The effective action U (P) for the Polyakov loop P(C) consists of a kinetic 
part and a potential part: 

U(P) = K(P)(dP) 2 + V(P). (18) 

In perturbation theory the kinetic term starts with the classical term (g 2 is 
the 't Hooft coupling): 

K(P) = ±(1+?K 1 (P) + ...), (19) 
7 



whereas the potential can be written as: 



V(P) = V 1 (P)+g 2 V 2 (P) + .... (20) 

Note the classical kinetic energy is proportional to 1/g 2 , compared to the 
(9(1) potential term coming from a fluctuation determinant.. This means 
the profile C(z) is slowly varying, ~ gz. 

We are therefore allowed to evaluate the effect of the fluctuations around 
a constant background value C of the phase of the Polyakov loop: 

A, = CS^o + gQji- (21) 

We will take Feynman background gauge. 

The first term in the potential is given by the fluctuation determinant: 



T 



V^PIC)) = — [-^log(-^ J D(C) 2 _)-rr 6 log(- J D(C) 2 _)-rr / log(i 7 . J D(C). 

(22) 

The ± sign means a commutator or anti-commutator. The sum over 
colour and momentum degrees of freedom is written as tr and is normalised 
by the three dimensional volume Vol. The suffix b (f) means doing the trace 
consistent with periodic (anti-periodic) boundary conditions. Tr sums also 
the spin degrees of freedom. 

The 7 matrices are taken care of by squaring the Dirac operator and 
correcting with a factor 1/2 in front of the logarithm. The trace over the 
Dirac indices gives a factor 4, so all in all: 

Tr f \og(i r D(C) + ) = 2tr f \og(-D(C) 2 + ). (23) 

The Z 2 symmetry leads to an effective potential with two degenerate 
minima, one where the Polyakov loop matrix takes the value 1 and one 
where it equals —1. In terms of the phase matrix C : 

C = 0, and C = 2ixT x ^diag(l, 1, 1, -1, -1, , -1) = 2irTY N/2 . (24) 

There are y entries with +1 and an equal number with —1. So P(C) = 
— 1 as promised. 

In both theories the tunnelling path is assumed to be the straight one 
connecting the two. That it describes a local minimum, i.e. a valley, is easy 
to establish. In gluodynamics the straight path hypothesis has been verified 
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for 3 and 4 colours to be indeed a global minimum 5 . For any number of 
colours in gluodynamics it leads to Casimir scaling, which is verified by lattice 
results [17]. 

Adding the two-index fermion deforms the potential. Only the Z 2 symme- 
try remains respected. We have checked that the straight line interpolation 
stays a local minimum. 

So along the tunnelling path we will take C to be of the form: 

C(q) = 2nTqY N/2 . (25) 

This is a straight path starting with P(C(q)) = for q = and ending 
in P(C(q)) = -1 for q = 1. 



3.2 Evaluation of the determinant 

The key element in eq.(22) is the covariant derivative. Since the back- 
ground is constant we expand in a plane wave basis. In this plane wave basis 
the trace is written like arf-l dimensional integral over momenta and a 
discrete sum over Matsubara frequencies: 

with 2e = 4 — d. Only from two-loop order on, the ji dependence will show 
up. 

To diagonalise the colour structure we use the Cartan basis with diagonal 
matrices A^, d = 1, TV — 1 and off-diagonal matrices A u with matrix 
elements (\ %1 )i m = ^6u6 jm . 

The off-diagonal basis elements are eigenvectors of any diagonal matrix 

C: 

[C,\ ii \± = {C i ±C j )\ ii . (27) 
In this basis the diagonalisation is trivial: 

D^C) ± Q^ =i(l, + (C l ±C ] )) (28) 

5 For the covering group of SO(2n + 4), having a Z4 centre group, the same hypothesis 
leads to inacceptable results 
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Because of the explicit form for C(q) (eq.(24) and (25)) the potential 
becomes: 



VM = 2 ) 2 [jj; log((Z„ + 2vrTg) 2 + ?)) - 2j£ \og(l + 2*T{q + l -)f + P) 

(29) 

The colour factor in front of the gauge contribution is 2(iV/2) 2 because 
the commutator in eq.(27) is only contributing if the indices i and j are 

in different sign sectors of Y N / 2 = l/2diag(l, 1, , 1, — 1, —1, , —1). The 

same colour factor results from the anti-commutator because in that case the 
non-vanishing contributions come from the same sign sectors. The factor 2 in 
front of the fermionic contribution comes from the spin, eq.(23), of our Dirac 
fermion. Since for fixed colour a Dirac fermion has twice as much degrees of 
freedom as the gauge boson this factor 2 is to be expected. This factor 2 is 
undone for each one of the irreducible representations, as we will see in the 
next subsection. 



3.3 Irreducible components of the Dirac fermion and 
equivalence 

Up till now we computed with a reducible two-index Dirac fermion. We 
now split the reducible representation ip^ into an anti-symmetric and a sym- 
metric representation with respectively |iV(iV ± 1) components. The colour 
factor 2(N/2) 2 in front of the reducible fermionic contribution in eq.((29) 
splits into two parts: |-/V(y ± 1) for symmetric and anti-symmetric repre- 
sentation 6 . 

For large- N both symmetric and antisymmetric representations contribute 
with the same colour factor (^) 2 . Combined with the spin factor 2 both for 
the symmetric and anti-symmetric representation we have from eq.(29): 

V 1± (q) = 2(^) 2 log((/ +27rTg) 2 +P))-(l±|)^ \og(l +2nT(q+ 1 -)) 2 +?)] 

(30) 

For large- N both the symmetric and antisymmetric representation con- 
tribute the same colour factor (y) 2 . The equivalence with the supersym- 
metric version is then very plausible since omitting the thermal boundary 



6 Note that the SU(2) antisymmetric representation does not contribute, as it is a colour 
singlet. 
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condition for the fermion, the boson and fermion contributions for V\ (q) can- 
cel out. 

Indeed our earlier result [6, 7] for the M = 1 evaluated for large-iV and 
the special channel k — y gives the same result for large- N. 

3.4 Potential and kinetic terms to two loop order 

First we discuss the potential. In general potentials as function of the 
Polyakov loop involve not only vacuum diagrams but also the renormalisa- 
tion of the loops. The latter gives rise to insertion diagrams and have been 
extensively discussed [18]. 

The two- loop diagrams for the potential are given in fig. (2). The only 
difference with the calculation in reference [7] is given by the fermionic loops 
in (b2) and (d). Since the latter come in linearly and the irreducible rep- 
resentations do not mix, we can calculate with the reducible representation 
and split the result later, as in the one loop case above. Moreover in diagram 
(d) only the derivative (with respect to the variable q) of the one loop result 
comes in. The inserted renormalisation of the Polyakov loop is to this order 
independent of the fermion content. So we know that the equivalence will 
work for diagram (d). 

As there are only four different sum integrals involved we list them here: 



The functions involved are proportional to the Bernoulli polynomials with 
the suffix indicating their order. They are shown in Appendix A. 

So the only diagram to analyse is diagram (b2). For k — y the result is 
for finite N: 




(31) 



(32) 



(33) 



(34) 
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Figure 2: One and two loop contributions to the effective potential. Continu- 
ous lines are bosons, broken lines are fermions and dotted lines are ghosts. In 
(a) the bosonic contributions, in (b) the fermionic contributions and (c) the 
ghost contribution are shown. In (d) the Polyakov loop (fat circle), renor- 
malized by one gluon exchange, is inserted into the sum of the loops (al), 
(bl) and (c), given by the shaded blob. 
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v - N2 

V2± — —r 



B 2 {qf + 2B 2 (0)B 2 (q) - (1 ± -){ - (B 2 (q + 1/2)) 2 + 2B 2 (q)B 2 (q + 1/2) + 
+ 2^B 2 (0)B 2 (q + 1/2) + B 2 (l/2)B 2 (q) - B 2 (l/2)B 2 (q + 1/2)) 



N 2 



B 2 (q + 1/2) 2 - 2B 2 (0)B 2 (q + 1/2)) } 



+ 



ABM (b-M -(l±^)B 3 (q+ 1/2)) 



The first two terms and the first term in the last line represent all the purely 
bosonic two-loop graphs in fig. 2, including the insertion diagram (d). The 
remaining terms in the curly brackets are the graph (b2) and the very last 
term is the contribution from the insertion graph (d). For N large we compare 
to the M = 1 result [7], using the notation D k (q) = B k (q) — B k (q + 1/2) to 
render the SUSY limit explicit: 



V 2 



(-) 2 
1 2 ; 



D2(q) + D 2 (0)) -D^Y + ABMDM 



(36) 



Indeed the equivalence is valid for the two loop potential. 

We could have seen this result almost directly by working in 't Hooft's 
double line notation. For the gluon we write: 



(Su'Sjj' - j-fSij'Si>j)/lii'- 



(37) 



with l 2 u , = (Z + Ci - cv) 2 + i 5 - 

Only in the two-index case the arrows on the fermion propagators run 
parallel. 

The colour shift is identical for adjoint or two-index case for colour charge 
Y N/2 , eq.(24). 

The 1/N part is the U(l) gluon we have to subtract from the first term. 
We can then sum over the indices as if they were independent as they would 
be in the U(N) case. 

In the Af = 1 case the adjoint fermion does not couple to the U(l) gluon. 
Hence, in the large- iV limit, the diagram (b2) in fig. (2), with the reducible 
two-index Dirac fermion and with the adjoint Dirac fermion running through 



(35) 
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(a) 



(b) 



Figure 3: Double line representation of the graph in fig. (2) (b2): equality 
in the large- N limit between an adjoint Dirac fermion (a) and a two-index 
Dirac fermion (b) running through the fermion loop. The difference in the 
(anti-)parallel arrows of the (adjoint) two-index fermions is irrelevant for the 
charge Fat/2 in the vertices. 

the loop are equal if the colour charge in the vertices is Yjv/2 as shown in 
fig. (3). The same equality holds when as in the J\f = 1 Major ana 

fermion is used, and in the two-index case an irreducible component appears. 



3.5 Kinetic term 

The remainder of this section is devoted to the one loop kinetic term 
shown in fig. (4). 

Clearly the first graph involves the renormalisation of the coupling. 
The beta function coefficients are 



b b = -11/3 
Po = b b + b f 



(38) 



In terms of the renormalized coupling g(T) of ref. [1]: 



1 



~9 2 {T) g 



9_ 

An 



- + ^(1/2)+ log 

£ 



TFT 2 



(39) 



we obtain for the kinetic term (ip(q) is the logarithmic derivative of the 
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Figure 4: The contributions to the one loop effective kinetic term. The 
wiggly lines are q lines. The shaded blob in (a) is the sum of one loop boson, 
fermion and ghost graphs. In panel (b) the renormalization of the Polyakov 
loop is shown. 



gamma function): 

Kl± - f?)Wfi-®y <w *(,) + «i-,))] + 



2 J ~g\T) \ V 4tt 

+ ((1±^)(6/W9 + ^)+^(-? + ^)) + y)+4)M- (4°) 

The last term is due to the insertion diagram (b) in fig (4). The ip function 
comes in through summation over the Matsubara modes: 

Eu7^ = r w?) ^ (1 - 9)) ' (41) 

n ' ' 

Let us compare this result for large N to the J\f = 1 result [7]: 



1 1 13 1 \ 

+ b f {^{q + -)+^ rea {-q+-)) + - + ^y 

Also the one loop kinetic term shows equivalence. Note: 



(42) 



the term ip reg {l/2 — q) in the Af = 1 case is defined by subtracting the 
pole part, 

^(1/2 - q) -+ Ae 9 (l/2 -q)= ^(1/2 - q) - 1/(1/2 - q). (43) 
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This pole originates in the n — — 1 Matsubara mode (see eq.(41) with 
argument q + 1/2 instead of q). Its subtraction is justified [7] because 
our fermion is a Majorana spinor and then of the two spin states with 
Matsubara frequency —1, one is not normalisable, and the other one is 
normalisable but does not contribute to the energy of the wall [19]. 

• The question now is how to regulate this pole in our two-index case, 
eq.(40). If the fermion is a Dirac fermion, with n — — 1, then there 
are four states. There are two finite bound state eigenspinors with zero 
energy. The other two are non-normalisable. This means that none of 
the n — — 1 fermion states contributes to the kinetic energy term, so 
that we can work with ip reg (l/2 — q) in eq.(40). 

• In the bosonic contribution there a is pole as well, at q — 0, but its origin 
lies in the thermal infra-red. It does not contribute to the tension, since 
it is folded with the square root of the lowest order potential B±{q). The 
latter is 0{q) near q — 0. 

We conclude that to two-loop order the effective actions of the two ir- 
reducible two- index fermions are identical to the Af = 1, up to terms of 
0(1/N). 

3.6 A note on perturbative equivalence for Af = 2,4 
SYM at high T 

At zero temperature the equivalence between Af = 2 and the two-index 
fermion theory becomes tenuous, because of the degeneracy of the ground 
state in the former. But at high temperature this degeneracy should vanish, 
as supersymmetry is now broken. 

So we felt justified to describe below the perturbative equivalence between 
them. In doing so we will also include Af = 4. 

Consider the action with two-index fermions in D = 6 or 10 dimensions 
and compare it to Af = 1 action in the same dimensionality. The two-index 
fermion has no constraint in 6 dimensions, whereas in 10 dimensions it obeys 
a Majorana (or alternatively a chiral) constraint. 

Upon dimensional reduction the Af — 1 D — 6 and D — 10 theories 
become Af = 2 and Af = 4 super Yang-Mills theories in D = 4. 

We consider these theories at high temperature and will assume that there 
is a deconfined groundstate with the scalars in the perturbative ground state. 
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Perturbation theory gives then identical results for the effective poten- 
tial U(P). The reason is that dimensional reduction does not change the 
integrals, only the counting of degrees of freedom. But the change in de- 
grees of freedom turns out to be the same in the two-index theory as in the 
corresponding SUSY theory. 

To one loop order the counting of the degrees of freedom due to the extra 
dimensions is simple. For the boson sector the comparison is immediate, 
for the fermion sector we get in 5 = 6 an eight dimensional spinor. Like in 
the four dimensional case (see eq.(30)) there are twice as many fermions as 
bosonic degrees of freedom, for a given colour. The same is true in ten dimen- 
sions, where the spinor is sixteen dimensional. This factor two is absorbed 
by going to the irreducible two-index fermion representations. 

For the one loop kinetic energy the same reasoning holds. And also for the 
insertion diagram (d) in fig. (2), because the renormalisation of the Polyakov 
loop stays unchanged in this order (Remember the Polyakov loop concerns 
only one polarization direction, the one in the periodic direction). 

In two loop order we need to analyse the diagram (b2) in fig. (2). All 
we need to do is convince ourselves that the equality in fig. (3) still holds 
in the 6 and 10 dimensional reductions of the respective theories. All what 
happens is that the fermion loop is corrected by not only a vector particle 
but also scalar particle exchanged. The only difference is factors of 1 ± ^ in 
the two-index case. 

We conclude that at the perturbative level there is, in addition to the 
Af = 1 case, equivalence between Af = 2,4 SYM and the corresponding 
"orientifold" theories. Note the equivalence is between the effective potentials, 
not only for its minimum, the tension, to which we turn now. 



3.7 Tension to two loop order 

As discussed in previous sections the tension equals the minimum of the 
effective action: 



NY AttT 



9 Jo 



dq 



V 2 - 



\?± ^~9 2 {^f + K 1± vl)+0(g 3 )] (44) 



All coefficients under the integral sign are those of eq.(30), eq.(35) and 
eq.(40), with the factor (^) 2 removed. 
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The resulting large- iV tension for the (anti-)symmetric two- index fermions 
is, from eq.(93) in [7]: 



a ±(T) 

{ N ?T2 



15 V m D (5) 



9 2 



(4tt)< 



{(-2.92683...) x b b (S) 



13. 



+ (3.27471...) xb f (S) + — } + Ax 5-2 
with A = (6 - 2)/2, m 2 D (5) = §g 2 T 2 and 



(45) 



b b (6) 
bf(5) 



11 5-4 
5-2 



(46) 
(47) 



As the ratio of Aj^r to T c is not (yet) known we cannot plot the tension 
eq.(44) as a function of T/T c . 



4 Orientifold planar equivalence for the wall 
tension - a nonperturbative proof 

Orientifold planar equivalence has been proved in different contexts and 
under different assumptions. A proof for local observables that assumes the 
perturbative expansion can be found in [9]. The first nonperturbative proof 
for the partition function and Wilson loops on the continuum was presented 
in detail in [10], in which the authors assumed that the fermionic determinant 
expansion in the worldline formalism is convergent (this assumption is exten- 
sively discussed in [20]). A proof for Wilson loops in the lattice-discretized 
theory, assuming the strong-coupling and large-mass expansions, was dis- 
cussed in [21]. Finally a nonperturbative proof for C-invariant operators was 
presented in [11], assuming that charge conjugation is not spontaneously bro- 
ken and that the coherent-state construction of [12] describes correctly the 
large- A limit of gauge theories. 

The setup we have in mind in this section is the one of the coherent states. 
This is particularly useful when dealing with expectation values of observ- 
ables expressed in the Hamiltonian formalism. A review of the coherent-state 
formalism and its application to orientifold planar equivalence is beyond the 
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scope of this paper. We will just summarize the main idea, and we will refer 
the reader to the literature [12, 11] for the details. 

The large- N limit of gauge theories is classical, in the sense that quantum 
fluctuations are suppressed. At every value of N it is possible to identify a 
particular overcomplete set of states (coherent states) in the Hilbert space. 
In the large-iV limit the set of coherent states becomes orthogonal and de- 
fines the classical phase space of the N = oo theory. Roughly speaking, the 
coherent states have an indetermination of order 1/N 2 , hence the indeter- 
mination vanishes in the large- N limit. This construction is analogous to 
the classical limit (h — > 0) for the coherent states of the harmonic oscilla- 
tor (indetermination of order K). Classical observables are operators in the 
Hilbert space that have the following properties: (a) their matrix elements 
with respect to the coherent states have a well defined large- iV limit, (b) they 
become diagonal in the basis of the coherent states in the large- iV limit. In 
the N = oo theory classical observables can be represented as functions on 
the classical phase space, and the commutator goes into the classical Poisson 
brackets. Roughly speaking, the classical observables are all the observables 
which are not sensitive to the quantum correlations between coherent states, 
and for which the large- iV factorization holds. Properly normalized Wilson 
loops (with possible insertions of electric fields and two-index fermions) and 
also the Hamiltonian are shown to be classical observables, while the vacuum 
is a coherent state. 

Consider now AdjQCD and OrientiQCD. The common sector is defined as 
the set of coherent states and classical observables that are charge-conjugation 
invariant. Assume that charge conjugation is not spontaneously broken in 
both theories, which means that the vacuum is in the common sector. In 
this setup, orientifold planar equivalence is stated as follows: 

1. a one-to-one map between the common sectors of AdjQCD and Orien- 
tiQCD exists; 

2. in the large- N limit the matrix elements of classical observables with re- 
spect to coherent states in the common sector are the same in AdjQCD 
and OrientiQCD; 

3. in particular the C-even spectrum of classical observables is the same 
in AdjQCD and OrientiQCD; 

4. since the vacuum is in the common sector, vacuum expectation values 
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of classical observables in the common sector are the same in AdjQCD 
and Orient iQ CD. 

As explained in Sect. 2, the domain wall tension is related to the thermal 
expectation value of the 't Hooft loop (k = N/2) in the deconfined phase (T 
is the temperature): 

V(L) = exp |^ J Tr(Y N/2 E)dS^ , (48) 

(V(L)) T = ^yTr V{L)e~ H ' T oc e~ L ^° , (49) 

where S is a rectangular L x x L y surface lying on the x — y plane at z — 0. 

One can be tempted to use the arguments above, in order to prove that the 
thermal expectation value of V(L) (and therefore the domain wall tension) 
is the same in AdjQCD and OrientiQCD in the large-iV limit. However, 
although the 't Hooft loop is charge-conjugation invariant, it is not a classical 
observable because it does not have a well defined large- N limit (it is e°( N •*). 
Moreover orient ifold planar equivalence cannot be straightforwardly exported 
to thermal expectation values. In fact in a thermal expectation value also 
C-odd states, which are not in the common sector, contribute. 

In order to bypass this problem, it is useful to express the 't Hooft 
loop thermal expectation value, as a zero-temperature magnetic correlator 
(Sect. 2) on a space with a compact dimension of length 1/T. We will 
prefer the lattice discretized formulation, since it avoids the technicality of 
dealing with discontinuous gauge transformations. We recall here the main 
formula (14): 

(V k (L)) = (0\eM-L y n t L )\0) , (50) 

where 1 0) is the vacuum of the Hamiltonian TL^. The Hamiltonians Hi and 
are defined respectively in eqs. (11) and (13). 

In the large L y limit, only the ground state \t) (in the C-even sector) of 
the twisted Hamiltonian % l L contributes to the numerator: 

(V n (L))~\(0\t)\ 2 e- L y E ^ , (51) 

where E t (L t ) is the energy of \t) with respect to the Hamiltonian % f L . If 
also L x is taken large, then the domain wall tension is recovered from the 
asymptotic behaviour of E t {L x ): 

E t (L x ) ~ aL x . (52) 
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The equality of the wall tension in the two theories will follow from the 
equality of E t (L x ). Consider the normalized twisted Hamiltonian: 

ht = ^ ■ (53) 
It can be written as the Hamiltonian in absence of twisting plus a correction: 

" twisted □ 

Separately, the normalized Hamiltonian l-i L /N 2 in absence of magnetic string 
and the correction (just a sum of normalized plaquettes) are classical observ- 
ables. Their matrix elements with respect to coherent states are finite in the 
large- N limit. Therefore h t is a classical observable itself. Moreover, since 
only the real part of the plaquettes appears, h t is also charge-conjugation 
invariant. Therefore it belongs to the common sector. By orientifold planar 
equivalence, its lowest eigenvalue E t (L x )/N 2 in the C-even sector is the same 
in AdjQCD and OrientiQCD. From Eq. (52): 

hm = ^ • (55) 

5 One-loop effective potential and the decay 
rate of the false vacua 

The one-loop effective potential for the Polyakov loop on thermal R 3 x Si 
and massive fermions in the antisymmetric representation is: 

2V3 J ^ ^ cos g(vi-Vj) (-l) g o-(gml3) cos g(v { + v j 



^'-S -EE T +M /EE 



1 ift 9>° % <3 9>0 



(56) 



x 2 



a(x) = -K 2 {x) . (57) 



Choosing the point Vi — ■ • ■ — Vjy — the energy density is (up to an 
additive constant that does not depend on k): 

_ V(2irk/N) 4N f N(N-l) ^ *(gmP)cos^ 

13 ' g>0 a 
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Since a(x) ~ -y^x 3 ^ 2 e x in the x — > oo limit, then for mf3 1 the only 
term contributing in the sum is g — 1 and: 

e(Jb) = v^ujm = _^mN^ {mP)3/ ^ C0S ^ (59) 

The energy difference between a metastable state and the stable state is: 

A £ (,) = eik) - £ (0) = ^f-^LzA^f,^ (, _ cos «) . 

(60) 

The large- N limit is taken in the following two cases: 



k is fixed: 



Ac(fc) = c(fc) - 6(0) = ^y|(m/3) 3 / 2 e-^A; 2 . (61) 

Ae(0) = iV 2 |§y|(m/3) 3 /V^sin 2 . (62) 

Consider now the following setup: the space is a R 2 x Si and antiperiodic 
boundary conditions have been chosen for the fermions along the compact 
direction. The system sits in a false vacuum, labeled by k, and then decays 
to the stable vacuum with k — 0. The case of SU (3) with sextet fermions 
was recently discussed in [22]. 

The decay rate of the k false vacuum is computed in terms of the action 
of the bounce, which is the solution of the classical Euclidean equations of 
motion, connecting the false vacuum to the true one in time. In the thin-wall 
approximation, the bounce looks like a bubble in the Euclidean spacetime. 
In four dimensions, it is a four-dimensional bubble. When a dimension is 
compactified with small size /3 = 1/T, it will look like a three-dimensional 
bubble of radius R. 

The action of the bubble is given by the sum of the action of the inner 
part of the bubble plus the contribution of the wall: 

S(R) = ~7rR 3 /3Ae k + AitR 2 a k (63) 
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where Ae k is the energy density difference between the false and true vacua, 
and <jfc is the domain wall tension. 

The maximum of the action is reached for: 

R = (64) 



S = T'wh (65) 



The thin-wall approximation is valid if Rfi ^> 1 where fi is the second 
derivative of the potential in the vacuum. This condition is fulfilled in the 
large- mass limit. 

The general formula for the decay rate of a false vacuum is: 

r ^^{ J T*wk] (66 > 

Consider the orientifold theory with a mass m for the fermions. When 
the mass becomes infinite, the fermions decouple, the theory becomes Zn 
symmetric and T = 0. When the mass m is very large but finite, the first 
nonzero contribution to T is obtained considering the first nonzero contribu- 
tion to Ae k (as computed in the section before): 

and the domain wall tension a k of pure Yang-Mills: 

n = k(N ~ k) wk- (68) 

The one-loop decay rate in this approximation (and every value of N) is: 

The large- N limit is taken in the following two cases: 
• k is fixed: 

T ^^<-T¥Njm^ e {mm ~ 3 > (70) 
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T k oc exp |- N ^-^ e W(mB)A (71) 

l fc ocexp| s . n4 ^ 34iV 2 ( 3 A) 3/2 e \ m P) J W 

6 Discussion 

In this paper we discussed theories with matter in the adjoint and two- 
index representations at high temperature. We focused on domain walls. 

The main results of our paper are: (i). a two-loop calculation of the 
tension of the domain wall the interpolate between the vacua with (P) = 1 
and (P) = —1, (ii). a comparison with the tension of the corresponding 
domain wall in a theory with adjoint fermions, (iii). a proof that at large- TV 
there is an exact equivalence between the domain walls of the two theories 
and (iv). a calculation of the decay rates of the false vacua in the "orientifold 
field theory". 

Our results suggest that planar equivalence holds not only in a given vac- 
uum, but also for objects that interpolate between distinct vacua. Moreover, 
planar equivalence holds not only for light (namely 0(1)) objects such as 
glueballs, but also for heavy (in the present case 0(N 2 )) objects such as 
domain walls. 

Concerning the vacuum structure of the large- N and large fermion mass 
orientifold theories: we learnt that those theories admit N — 2 false vacua 
(that become true vacua when the fermion mass becomes infinite and de- 
couples). Those vacua have a very narrow width in the large- N limit: they 
decay exponentially either as exp — iV 3 (70) or exp — N 2 (71), depending on 
the way that the large- N limit is taken. 

There are several interesting future directions to explore. It will be inter- 
esting to study other theories with a Z 2 centre, such as theories based on an 
orthogonal gauge group and in particular to study their gravity dual. Such 
theories contain unoriented strings. 

It is also interesting to study, from the string dual side the existence and 
decay of the false vacua into the true vacua of the theory. 

Finally, as so little is known about deconfining domain walls from lattice 
simulation, it will be interesting to carry out simulations of domain walls of 
theories with fermions in higher representations on the lattice. 
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A Details about the two loop calculation 

The Bernoulli polynomials [23] are related to the sum integrals in eq. (34) 

as: 

B d=4 (x) = \ 2 T 4 B 4 (x) (72) 
B 3 (x) = \nT 3 B 3 {x) (73) 
B 2 {x) = \t 2 B 2 {x) (74) 
B x {x) = ~B x {x) (75) 

and finally (e(x) is the sign function): 

B 4 (x) = x 2 (l - \x\) 2 (76) 
3 1 

B 3 (x) = x 3 - -e(x)x 2 + -x (77) 

B 2 (x) = x 2 -e{x)x+ l - (78) 

B 1 (x) = x-±e(x) (79) 

With this definition the relations (75) are valid in the interval — 1 < x < 1. 

Even (odd) polynomials are even (odd) in x. From the original definition 
(34) all the B(x) are periodic mod 1. 

We now turn to an analysis of the two-loop formulae for the potential 
V 2 ±. 

If the weights are all the same, then the total multiplicity is N 3 by adding 
up the entries in the multiplicity column. 

We took the indices i,j,m to be independent. Hence we have to correct 
by subtracting the diagram with the U(l) gluon which has multiplicity N so 
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Table 1: Counting of the multiplicity of the two-loop planar diagram, [k] is 
the Y k sector with k entries k — N, and [N — k] is the sector with N — k entries 
k. i,j and m are the indices of the index cycles discussed in the text. The 
function B in column 4 equals B?,{x) = X i ((^ + 2nTx) 2 + l 2 )^ 1 , as appears 
in eq.(34). The full result for the diagram b2 results from multiplying the 
multiplicities and the weights for every row. 



[k] 


[N-k] 


multiplicity 


weight 


mj 

i 


i 

mj 


k' 2 (N — k) 
k(N-k) 2 


B(2(r - 1)9 + i)S(i) - B(q)(B(2(r - l)q + ±) + fl(±)) 
5(2rg + ±)fl(±) - B(q)(B(2rq + \) + B(±)) 


j 

mi 


mi 
j 


k(N-k) 2 
k 2 (N -k) 


B(2rq + i)S(i) - B(q)(B(2rq + ±) + 5(A)) 
S(2(r - l)g + |)5(|) - B(g)(B(2(r - 1) + §) + B{\)) 


m 
U 


U 
m 


k(N-k) 2 
k 2 (N - k) 


B(2r - Iq) + A) 2 - 25(0)5(2r - Iq) + ±) 
5(2r - l)q + \) 2 - 25(0)5(2r - l)q + \) 


ijm 


ijm 


k' A 

(N - kf 


B(2(r - l)q + |) 2 - 2B(0)B(2(r - l)q + ±) 
5(2rg+±) 2 -25(0)5(2rg+±) 



the total multiplicity at g = is N(N 2 — 1). The weight of the £7(1) gluon 
contribution combined with the multiplicity is: 



1 

iv 2 



+ 



k 2 ((B(2(r - l)q + i) 2 - 2£(0)5(2(r - l)q + 1)) 

(iV - k) 2 ({B{2rq + i) 2 - 2B(0)B(2rq + I)) 

+ 2k(N - k) (B((2r - l)q + i) 2 - 2B(0)B((2r - l)q + i)))] . (80) 

The factor ^ comes from the U(l) part of the double line gluon propa- 
gator: 

Its contribution is 0(1/ N 2 ) smaller than the leading term, because of the 
explicit 1/N factor in the propagator and one index loop less in the diagram. 
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